Mesoscopic Correlation with Polarization of Electromagnetic Waves 
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Mesoscopic correlations are observed in the polarization of microwave radiation transmitted 
through a random waveguide. These measurements, supported by diagrammatic theory, permit 
the unambiguous identification of short, long, and infinite range components in the intensity corre- 
lation function, as well as an additional frequency-independent component. 
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Mesoscopic transport of both classical and quantum 
mechanical waves is characterized by the degree of non- 
local intensity or current correlation which reflects the 
closeness to the Anderson localization threshold 0, 
Correlation is generally treated using scalar waves with 
the electromagnetic (EM) polarization or electron spin 
accounted for by doubling the density of states. The 
vector nature of EM radiation is exhibited in coherent 
backscattering the photon Hall effect^3, and the 
memory of rotation of incident polarization ja| . These po- 
larization effects reflect short-range fleld correlation and 
can be obtained within the fleld factorization approxi- 
mation applicable to Gaussian statistics, as opposed to 
mesoscopic intensity correlation. 

In this Letter, we present measurements and calcu- 
lations of intensity correlation with polarization rotation 
that go beyond the fleld factorization approximation. Po- 
larization is thereby shown to be an essential mesoscopic 
variable complementary to space, frequency, and time. 
Unlike these variables, however, the polarization vari- 
able is of flnite range with only two independent values. 
Consequently, the fleld correlation function with polar- 
ization rotation has a simple form independent of scatter- 
ing strength. It is determined exclusively by the cosine 
of the angular shifts in the polarization of the electric 
fleld at the source and detector and vanishes when either 
polarization is shifted by 90°. This enables an unam- 
biguous separation of the intensity correlation function, 
CJnto short (Ci) 0,0, long (C2) and inflnite (C3) 

|lfl IllL range components, which determine fluctu- 
ations in intensity, total transmission, and conductance, 
respectively. The frequency correlation function of each 
of these universal components can thereby be found. In 
addition, the presence of a frequency-independent com- 
ponent is identifled. This term may correspond to a 
nonuniversal "Cq" contribution, which reflects the scat- 
tering environment at the input and output surfaces . 

We begin with the conjecture, borne out by diagram- 
matic calculations and measurements reported below, 
that the structure of the correlation function of normal- 
ized intensity with regard to rotations in the polariza- 
tion of source and detector, C {/S.6 $ , o) , is analogous 



to that of the correlation function with displacement of 
the source and detector, C(Ar5', Aru) In each com- 
ponent, C depends upon a given degree of freedom only 
through its dependence upon the square of the fleld cor- 
relation function, F = \Fe\^ ■ When expressed in terms 
of polarization rotations of the source and detector, this 
gives, 

C{Ms.Md) - F(Ms)F{Md) 

+ A'^[F{Ms)+F{Md)] 

^F{Ms) + F{Md) 
F{Ms)F{Md)]. (1) 
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The leading contributions to A'2 and A'^ are of order 1/g 
and 1/5^, respectively, where g is the dimensionless con- 
ductance. 

The field correlation function with polarization rota- 
tion can readily be found for multiply scattered waves 
when the transmitted fleld is completely depolarized. In 
this case, the average intensity is independent of the po- 
larization of the fleld at the source or detector and cross- 
polarized flelds are uncorrelated. When the polarization 
of the detector is rotated by A0£) from an initial direction 
the detected fleld E(0£) -I- A0£)) may be expressed in 
terms of the vector sum of the field along 9^^ E(0£)), and 
the field perpendicular to this direction, E(6'£) -I- 90"), 
as follows, E(0D + Md) = E(6'd) cos A6'£, + ^{Od + 
9O°)sinA0£). The field correlation function for the nor- 
malized average intensity, {\Ej{6o)\'^) = 1, is therefore, 
Fe{Md) = {'E{eD)E*{9D + A9d)) - cosA6'i5. Here 
(...) denotes the average over an ensemble of realizations. 
A similar argument gives the field correlation function of 
the detected field when the source polarization is rotated 
by A9s, Fe(A6's) = cos A6'5. The total field correlation 
function with polarization rotation is therefore given by 
FEiAOs^Aeo) = FE{A0s)FEiA0D) - cosA^s cosAfl^^. 

The contributions to C in Eq. (1) can be expressed in 
terms of the product or sum of F{A9s) and F{A9d) or 
of a constant. Since these three types of terms are asso- 
ciated, respectively, with short, long, and infinite range 
correlation, it is convenient to group these terms sepa- 
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FIG. 1: Schematics of the experimental setup: S (source) and 
D (detector) microwave horns are positioned symmetrically in 
front and behind the sample and can be rotated about their 
axes (z-axis) thus rotating the polarizations of the incident 
and detected waves. 



rately. We therefore write, 

C{Aes,AeD) = [1 + Ai)cos^Aescos'^AeD 

+ A2{cos^Aes + cos^AeD) +A3.{2) 



For the case, in which nonuniversal "Co" contributions to 
C are negligible, Ai = A'^, A2 ^ A!^ + A!^, and A3 = 
The analysis of Ref. has been generalized to vector 
waves in a detailed diagrammatic study 15]. The rela- 
tion A\ = A3, which follows from Eq. (1), can be under- 
stood from the topological structure of the underlying di- 
agrams, and the equal weights of the two Ai terms follow 
from reciprocity. In the diffusive regime in the absence 
of absorption, to order 1/(7^, 



(3) 



where g — 2 y. Ak'^i/SnL, which includes both polariza- 
tions. Here, k is the wave number, £ is the transport 
mean free path, A and L are the cross-section area and 
length of the sample, respectively. Values of A2 and A3 
slightly depend on absorption [ij, 113 • In our exper- 
iment, L/La = 3.6, where La is the absorption length. 
This gives A'^ = 0.57/g and A^j = O.llS/g^. 

Measurements of the dependence of the field and in- 
tensity correlation functions upon polarization rotation 
for microwave radiation transmitted through a random 
dielectric sample are made with identical conical horns 
positioned 40 cm in front and behind the sample (Fig. 1). 
Linearly polarized microwave radiation is launched from 
horn S and a linearly polarized component of the trans- 
mitted field is captured by horn D. The polarization of in- 
cident and detected waves can be rotated in the ccy-plane 
by rotating the corresponding horn about its axis {z- 
axis). Because of the conical shape of the horns, the an- 
gular distribution of the incident and detected radiation 
do not depend strongly on the rotation of the horns. As 
a result, the variation in the detected field at a given fre- 
quency for a given sample realization is mainly due to the 
shift in the polarization of the field at the source and de- 
tector. The sample is composed of 0.95-cm-diameter alu- 
mina spheres with refractive index 3.14 embedded within 
Styrofoam spheres of diameter 1.9 cm and refractive in- 
dex 1.04. The sample with an alumina volume fraction 
of 0.068 is contained in a 7.3-cm-diameter copper tube 



of length L = 90 cm. Measurements for random ensem- 
bles are obtained by momentarily rotating the sample 
tube about its axis to create new sample realizations, af- 
ter which spectra are taken for each orientation of the 
polarizers. Measurements are carried out in an ensem- 
ble of 12,000 sample realizations for 2 orientations of the 
source polarizer rotated by 90° and 7 orientations of the 
detector polarizer rotated in steps of 15*^, and also in 
an ensemble of 32,000 realizations for 2 orientations of 
each polarizer rotated by 90". The transmitted field is 
measured in steps of 1 MHz from 16.95 to 17.05 GHz 
using a Hewlett-Packard 8772C vector network analyzer. 
This frequency range is centered at the peak of the fourth 
Mie resonance of the alumina spheres |l6| and is much 
narrower than the width of the resonance so that propa- 
gation parameters within this range are nearly constant. 
At z/ = 17 GHz, the number of transverse channels in 
the sample tube is = Afc^/27r = 84. The transport 
mean free path estimated from Mie theory is £ = 2.34 
cm, giving g « 4:N£/3L = 2.91. A fit of diffusion theory 
to the measured field correlation function with frequency 
shift gives D = 8.0 cm^/ns for the diffusion coefficient 
and La = 24.9 cm for the absorption length. 

The real and imaginary parts of -^^(O, AOd) obtained 
by averaging over frequency for an ensemble of 12,000 
sample configurations are displayed by the crosses in 
Figs. 2a and 2b, respectively. The dashed curves rep- 
resent the predictions — cosA^d for the real part and 
for the imaginary part. The orientation of the horn 
polarizations for the initial and final data points are in- 
dicated by the arrows. In addition, we find Fe(90, AOd) 
to vanish for all the values of A9d- Though the measure- 
ments are consistent with the predicted behavior of Fe, 
there is an increasing systematic deviation from theory 
with increasing Ads or A9d. This deviation arises be- 
cause, in addition to rotation of the polarization, there is 
a residual variation of the spatial intensity distribution 
on the sample, as the horn is rotated. 

The intensity correlation functions C{0,AdD) and 
C(90, A9d) for the same data set are shown by the blue 
crosses in Figs. 3a and 3b, respectively. In addition, 
F{A9d) is shown in Fig. 3a by the red crosses, which 
lie close to the solid curve of cos^Ado- A Legendre 
polynomial fit of C(0, A9d) and C(90, A9d) shows that, 
with high accuracy, they are of the form a + 6cos^A6'd, 
where a and h are constants, in agreement with Eq. (2). 
The dashed curve in Fig. 3a is a fit to the first 4 data 
points of C(0, A9d) of the function a + hcos^A9m where 
a = A2 -I- A3 and h = 1 -I- Ai -I- A2. Only the first four 
points are used because the systematic deviation in C due 
to the redistribution of intensity with rotation cannot be 
ignored for larger polarization rotations. Though uncer- 
tainty is introduced into Ci and C2 by the horn rotation, 
there is no such effect on C3 because it is independent of 
both polarization and displacement. Thus the measured 
value for A3 is not affected by the horn rotation. It is 
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FIG. 2: (a) Real and (b) imaginary parts of the field correla- 
tion function Fe{0, A6d) plotted versus the angular shift of 
the polarization of the detected field, A9o, as the polarization 
of the incident wave is kept unchanged, AOs = 0. The arrows 
indicate orientations of the source (S) and detector (D) po- 
larizers for the initial and final data points. The dashed lines 
represent the theoretical prediction, Fe{0, A6d) = cos(ASo). 



found directly from C(90,90) for an ensemble of 44,000 
realizations, ^3 = 0.029 ±0.001. From the fit, we obtain 
Ai+A2^ 0.286 ± 0.005 and A2 + ^ 0.293 ± 0.004. 
Hence A2 = 0.264 ± 0.006 and Ai = 0.022 ± 0.010. The 
dashed curve in Fig. 3b is As + A2 cos^A^d with the A2 
and A3 found above. The difference between the exper- 
imental data points and the curve represents the loss of 
correlation due to the redistribution of intensity with the 
horn rotation. 

The coefficients Ai for arbitrary frequency shift can 
be determined from spectral measurements for three ori- 
entations of the source and detector as follows, A3 = 
C(90,90), A2 = C(0,90) - A3, and Ai = C(0,0) - 
|Fe(0,0)P - 2A2 - A3. We have utilized the above 
procedure in the ensemble of 44,000 sample realiza- 
tions, to find Ai{/S.v) for frequency shifts /S.v up to 100 
MHz, that is almost 100 times the Thouless frequency 
Ai^T/i = QD/2ttL'^ sa 1 MHz. To compensate for the 
loss of correlation in C(0, 90) due to the redistribution of 
intensity, the corresponding frequency correlation func- 
tion was multiplied by a constant to make it equal to 
the A2 -I- A3 at zero frequency shift. The resulting de- 
pendencies are shown by the squares in Fig. 4. A2(Ai^) 
and A'i{Av) are predicted 0| to fall asymptotically as 
A2(Az/) ~ l/(A:/)i/2 and As{l^v) - l/(Ai/)3/2. We 
find, however, that A3 exhibits a nonzero correlation for 
Ai/ 3> Az^Th, that largely exceeds error bars, and that A2 
is shifted from the predicted curve by a nearly constant 
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FIG. 3: (a) Intensity correlation functions C and F plotted 
versus AQd for Ads = 0. The solid line is the prediction, 
F{ABo) =cos'^{A0d). The dashed line is the fit of C(0,Ae'i3) 
to the first 4 data points, (b) C versus Ado for AOs = 90*^. 
The dashed line is the prediction A3 + A2 cos^A^d with the 
A2 and A3 found from the fit in (a). 



value. Figs. 4a and 4b show that agreement with theory is 
obtained for A2,3, if frequency-independent contributions 
with the values A2(oo) — —0.015 and A3 (00) — 0.007 are 
incorporated, respectively, in A2 and A3. We do not ob- 
serve an asymptotic background in Ai (Fig. 4c). The 
observed Ai(Ai/), however, falls faster than predicted 
by diagrammatic theory. To calculate Ai(A/), we have 
modified the theory for A3 (A/) 0. 

After subtracting the asymptotic values we find 
A3(0) - A3(cx3) = A^(0) = 0.022 ± 0.001 and A2(0) - 
A2(oo) = A^(0) -I- A(j(0) = 0.279 ± 0.004 at zero fre- 
quency shift. Note that A3(0) = Ai(0) holds within 
the error bars. Using Eq. (3) the value for A3(0) yields 
g — 2.27 ± 0.05 for the dimensionless conductance. This 
differs from the estimate, reflecting effects of localiza- 
tion, internal reflection, and finite scatterer density. Us- 
ing this value, Eq. (3) predicts A2(0) = 0.252, and hence 
A2(0) -I- A3(0) = 0.274, which is also consistent with ob- 
servations. 

It is tempting to associate the frequency-independent 
background term with the Cq term of Ref. A pre- 

liminary calculation for point dipoles in a tube geom- 
etry with N transverse channels shows Cq to be 
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FIG. 4: Coefficients A3 (a), A2 (b), and Ai (c) piotted versus 
frequency sfiift. Tlie iong-dashed curves represent predictions 
of mesoscopic tlieory [17] for A3 (a), A'2 (b), and Ai (c) with 
g = 2.27 and L/La = 3.6. Tlie frequency-independent contri- 
butions witfi tlie values A3 (00) = 0.007 and A2{oo) — —0.015 
have been added, respectively, to the predicted A3 and A'2 to 
obtain the solid curves and agreement with the data. 



independent of frequency shift, positive and of order 
— 0.012, and with contributions to A2 and ^3 
and not to Ai. This is consistent with the observed 
value A3 (00) = 0.007, but not with the negative value 
A2{oo) — —0.015. A vector Mie theory for Co along with 
measurements in samples with different longitudinal and 
transverse dimensions may shed light on this issue. 

In conclusion, mesoscopic correlations have been found 
in the polarization of microwave radiation transmitted 
through a random waveguide. The intensity correlation 
function with polarization rotation is shown to be com- 
posed of the product and sum of the square of the field 
correlation function and of a constant term. These three 
types of terms are associated, respectively, with short, 
long, and infinite range correlation. The variation with 
frequency shift of each of these terms has been measured. 
The long and infinite range terms are found to include 



a frequency-independent component, which may be re- 
lated to the short-duration interaction of the wave with 
regions of limited spatial extent near the input and out- 
put surfaces of the sample. The simple form of the three 
contributions to the intensity correlation function with 
polarization rotation allows an unambiguous separation 
of correlation into short, long, and infinite range con- 
tributions. This method is readily applicable to optical 
measurements and can be used to determine the degree of 
intensity correlation in a sample and hence the closeness 
to the localization threshold. Analogous correlations can 
also be expected in electron spin propagating in a random 
potential. 
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